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Abstract. This contribution presents a review of our present theoretical as well as experimental knowledge 
of different fluctuation observables relevant to nuclear multifragmentation. The possible connection between 
the presence of a fluctuation peak and the occurrence of a phase transition or a critical phenomenon is 
critically analyzed. Many different phenomena can lead both to the creation and to the suppression of 
a fluctuation peak. In particular, the role of constraints due to conservation laws and to data sorting is 
shown to be essential. From the experimental point of view, a comparison of the available fragmentation 
data reveals that there is a good agreement between different data sets of basic fluctuation observables, if 
the fragmenting source is of comparable size. This compatibility suggests that the fragmentation process 
is largely independent of the reaction mechanism (central versus peripheral collisions, symmetric versus 
asymmetric systems, light ions versus heavy ion induced reactions). Configurational energy fluctuations, 
that may give important information on the heat capacity of the fragmenting system at the freeze out 
stage, are not fully compatible among different data sets and require further analysis to properly account 
for Coulomb effects and secondary decays. Some basic theoretical questions, concerning the interplay 
between the dynamics of the collision and the fragmentation process, and the cluster definition in dense 
and hot media, are still open and are addressed at the end of the paper. A comparison with realistic models 
and/or a quantitative analysis of the fluctuation properties will be needed to clarify in the next future the 
nature of the transition observed from compound nucleus evaporation to multi-fragment production. 

PACS. 24.10.Pa - 24.60. Ky - 25.70.Pq - 64.60.Fr - 68.35.Rh 
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1 Fluctuations and phase transitions 



Since the first inclusive heavy ion experiments, multifrag- 
mentation has been tentatively associated to a phase tran- 
sition or a critical phenomenon. This expectation was trig- 
gered by the first pioneering theoretical studies of the nu- 
clear phase diagram which contains a coexistence re- 
gion delimited, at each temperature below an upper criti- 
cal value, by two critical points at different asymmetries [21 

Even more important, the first exclusive multifragmen- 
tation studies have shown that multifragmentation is a 
threshold process occurring at a relatively well defined 
deposited energy [4, 5, 6, 7 . The wide variation of possible 
fragment partitions naturally leads to important fluctua- 
tions of the associated partition sizes and energies. 

Different observables have been proposed to measure 
such fluctuations. Using the general definition of the n — th 
moment as 



the variance of the charge distribution is measured by the 
second moment M2 or by the normalized quantity [8]: 
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The root mean square fluctuation per particle 
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of the distribution of the largest fragment Z,„ detected in 
each event completes the information. We will also con- 
sider the total fluctuation 



y2 



and the fluctuation 
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of the configurational energy per particle associated to 
each fragment partition (fc) 
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where is the multipUcity of event fc, BE is the ground 
state binding energy of each fragment, and — is 
the average interfragment distance at the formation time. 
The quantities Aq,Zo in eqs.([3]),([5|) represent the recon- 
structed charge and mass of the fragmenting system, Zq ~ 

In a simple statistical picture the fluctuation of any 
observable can be related to the associated generalized 
susceptibility by 



(7) 



where A is the intensive variable associated to the generic 
observable A. Since the intensive variable associated to a 
particle density N/V is the susceptibility x = d{N)/d^, 
then the large variance of the charge distribution observed 
in multifragmentation experiments could be connected to 
the diverging critical point fluctuation which would signal 
a diverging susceptibility and a diverging density correla- 
tion length. The apparent self-similar behavior and scaling 
properties of fragment yields^ tends to support this in- 
tuitive picture. 




Fig. 1. Second moment of the size distribution (see [TT] for a 
precise definition) as a function of the average cluster multi- 
plicity for the three-dimensional percolation model for different 
lattice sizes. Figure is taken from 8 . 



1.2 Thermal invariance properties 



1.1 Finite size effects 



Many different effects can however blur this simple con- 
nection. First of all, since fragmenting sources cannot ex- 
ceed a few hundred nucleons, we have certainly to expect 
finite size rounding effects, which smooth the fluctuation 
signal [5|. Not only the transition point is expected to be 
loosely defined and shifted in the finite system as shown 
in the three-dimensional percolation model in Figure [H 
but also the signal is qualitatively the same for a critical 
point, a first order transition or even a continuous change 
or cross over. 

Finite size effects have other consequences on the dis- 
tribution than the simple smoothing of the transition. 
It has been shown on different model calculations that 
the presence of conservation constraints as well as the 
use of different event sorting procedures can sensibly dis- 
tort the fluctuation observables. To give a simple exam- 
ple, the presence of a peak in the largest fragment's size 
fluctuation as a function of the energy deposit is trivially 
produced by the baryon number conservation constraint 
which forces this fluctuation to decrease with increasing 
average multiplicity [10]. In the case of a genuine critical 
behavior as for the percolation model, the fact of sorting 
events according to the percolation parameter p or accord- 
ing to some other correlated observable, as for instance 
the total cluster multiplicity, modifies [lOlfS] the behavior 
of 7712, 72, and all other related moments [11 j measuring the 
fluctuation properties of the system. All these effects can 
be understood in the general framework of the non- equiv- 
alence of statistical ensembles for finite systems, which we 
will discuss in the next chapter. 



Another problem when trying to connect a fluctuation 
peak to a phase transition or a critical behavior in a fi- 
nite system is given by the possible existence of thermo- 
dynamic ambiguities. It has been observed by different 
independent works that in the framework of equilibrium 
fragmentation models the fluctuation behavior is qualita- 
tively independent of the break up density [T^13[[H1[T5] . 
An example is given in Figure [2l which gives the second 
moment of the charge {S2 = M2 — Mf) and of the en- 
ergy (Cy) distribution as a function of temperature in the 
Lattice Gas model for different break-up densities in the 
subcritical regime. 

A peak in the fluctuation observables can be seen at all 
densities, at a temperature which is systematically below 
the critical temperature of the system and close to the first 
order transition temperature in the thermodynamic limit. 
A similar behavior has been observed in different fluctu- 
ation observables and also at supercritical densities along 
the Kertesz percolation line, where the system does not 
present any phase transition. Table [1] gives, as a function 
of the lattice size, the inverse temperature at which the 
variable S2 shows a maximum in the three-dimensional 
IMFM model[T3] at different densities. As a general state- 
ment, the fluctuation peak as well as the global scaling 
properties of the size distribution [ISlfTB] in these models 
can be found along a curve in the T(p) diagram passing 
through the thermodynamic critical point but extending 
in the subcritical as well as supercritical region[l3- The 
subcritical behavior can be understood as a finite size ef- 
fect, when the correlation length, close to the first order 
transition point, becomes comparable to the linear size of 
the system, while the supercritical behavior is linked to the 
definition of clusters in dense and hot media[T7]. For the 
subcritical region, a clusterization algorithm has been sug- 
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Fig. 2. Second moment of the charge {S2) and of the energy 
(Cv) distribution as a function of temperature in the Lattice 
Gas model for different densities for a system of linear dimen- 
sion L=7. Arrows: first order transition temperature in the 
thermodynamic limit. Figure taken from ref. [14j. 



Table 1. Inverse temperature at which the second moment 
S2 = M2 — Ml is maximal for different densities and lat- 
tice sizes in the three-dimensional IMFM model. Taken from 
ref. [13] 
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gested to eliminate such behaviors in Ising simulations [9]. 
The possible pertinence of all these observations to exper- 
imental data is still a subject of debate, and essentially 
depends on the relationship between the measured clus- 
ters and the cluster definitions of the models. 

Last but not least, the presence of different timescales 
in the reaction |T8 l ll9 1 [20] and the dynamics of the frag- 



mentation process may have important effects in the quan- 
titative value of charge partition fluctuations [2 Ij . as we 
will discuss in the last chapter. 

For all these reasons, it is clear that the well docu- 
mented presence of a fluctuation peak in the measured 
charge distributions [7 cannot be taken as such as a proof 
of a critical behavior and/or phase transition. In order 
to connect the fluctuation behaviour to a phase transition 
and to conclude on its order it is indispensable to compare 
with models and/or to quantify the fluctuation peak. 



2 Theory 

2.1 Fluctuations and constraints 

It is clear that fluctuations on a given observable A will 
be suppressed if a constraint is applied to a variable cor- 
related to A. This trivial fact has a deep thermodynamic 
meaning and is linked to the non equivalence of statistical 
ensembles in finite systems [22] . Indeed the basic statistical 
relation between a fluctuation and the associated suscep- 
tibility eq.© is only valid in the ensemble in which the 
fluctuations of A are such as to maximise the total entropy, 
under the constraint of (A) ("canonical" ensemble). The 
thermodynamics in the ensemble where the generic ob- 
servable A is controlled event by event (" microcanonical" 
ensemble), or in the ensemble where a a is externally fixed 
("gaussian" ensemble [23]) is a perfectly defined statisti- 
cal problem, but the thermodynamic relationships have 
to be explicitly worked out [24]. As an example we show 
in Figure [3] the correlation between the size of the largest 
cluster Abig and the total energy in the isobar Lattice Gas 
model[24] at the transition temperature. The presence of 
two energy solutions at the same temperature and pres- 
sure clearly shows that the transition is in this case first 
order pS]. The Abig fluctuation properties are very differ- 
ent in the canonical ensemble (left part) and in the mi- 
crocanonical ensemble (right part) at the same (average) 
total energy. Because of the important correlation between 
the total energy and the fragmentation partition, fragment 
size fluctuations can be compared only for samples with 
comparable widths of the energy distribution. 

From the experimental viewpoint, different constraints 
apply to fragmentation data and have to be taken into 
account. Apart from the sorting conditions [7] , the coUi- 
sional dynamics can also give important constraints to 
the fragmentation pattern (e.g. flows, deformation in r- 
space and p-space). This means that fluctuations have to 
be compared with calculations performed in the statisti- 
cal ensemble corresponding to the pertinent experimental 
constraints [26, . 



2.2 Fluctuations and susceptibilities 

In the last subsection we have stated that a connection be- 
tween a fluctuation and the associated susceptibility can 
always be in principle worked out if the constraints acting 
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Fig. 3. Center: correlation between the largest fragment's size 
and the total energy in the isobar Lattice Gas model close to 
the transition temperature, for a system of 216 particles. Left 
side: projection over the Aug direction. Right side: same as 
the left side, but only events within a narrow energy interval 
around the average energy have been retained. 



on the observable are known. In the case of sharp con- 
straints {e.g. fixed total mass, charge, deposited energy), 
the connection between the fluctuation on a variable cor- 
related to the constraint {e.g. size or charge of the largest 
fragment, configurational energy) and the associated sus- 
ceptibility are in many cases analytical 27,28,29J. If a con- 
servation constraint A — Ai + A2 = est applies and the 
system can be splitted in two statistically independent 
components such that W{A) = W{Ai)W{A2), then the 
partial fluctuations are linked to the total susceptibility 

by 



X 



'ref 



(8) 



where = d\.Wi, cr^g^ is the fluctuation of Ai in the 
ensemble where only the average value {A) is constrained, 
and we have approximated the distribution of Ai with a 
gaussian[24]. The case of the total energy constraint has 
been particularly studied in the literature. Indeed the total 
energy deposit can be (approximately [3D]) measured event 
by event in An experiments, allowing to experimentally 
construct a microcanonical ensemble by sorting. For clas- 
sical systems with momentum independent interactions 
the potential energy fluctuation aj at a fixed total energy 
is linked to the total microcanonical heat capacity by 



1 



(9) 



where Ck,C are the kinetic and total heat capacity, cr^ = 
and tT^^„ = c^T^ is the kinetic energy fiuctuation in 
the canonical ensemble. Apart from the microstate equi- 
probability inherent to all statistical calculations, the above 
formula is obtained in the saddle point approximation for 
the partial energy distributions. The contribution of non 
gaussian tails can be also analytically worked out [28] and 
has been found to be negligible in all theoretical as well as 
experimental data samples analyzed so far|24j. An exem- 
ple of the quality of the approximation is given in Figure |4] 
which gives the temperature, normalized potential energy 
fiuctuation and heat capacity in the isobar Lattice Gas 
model for a system of 108 particles. 
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Fig. 4. Temperature, normalized binding energy fluctuation 
and heat capacity in the microcanonical isobar Lattice Gas 
model as a function of the total energy for a system of 108 
particles. In the lower panel the heat capacity estimated from 
fluctuations via eq.Q (dots) is compared to the exact expres- 
sion from the entropy curvature (line). Figure is taken from 

3 Experiment 

3.1 Effect of the sorting variable 

In this section we turn to compare different sets of experi- 
mental data available in the literature. A special attention 
has been paid by different collaborations to the largest 
fragment fluctuation am eq.Q and to the 72 observable 
eg . (|^ ^5> ,32..33.,34.,g5ll5S] . For all data sets of comparable 
total size these observables, as well as the others we will 
show in the next subsections, show a well deflned peak 
at comparable values of the chosen sorting variable. This 
is an important and non trivial result considering that 
data are taken with different apparata and the multifrag- 
menting systems are obtained with very different reaction 
mechanisms. The effect of the sorting variable is explored 
in table [21 that gives the maximum value of 72 and am 
with different data sets sorted in bins of the total mea- 
sured bound charge Ztound, total measured charged parti- 
cles multiplicity m, or calorimetric excitation energy [5I?j. 
Even if the systematics should certainly be completed and 
errors should definitely be evaluated, we can observe from 
table [2] that different data sets show a reasonable agree- 
ment when the same sorting is employed. 

We can also note that a higher 72 is systematically 
obtained when data are analyzed in bins of total charge 
multiplicity, with respect to a sorting in Zbouiid- This can 
be qualitatively understood if we recall that 72 measures 
the variance of the charge bound in fragments, and this 
quantity is obviously strongly correlated with Zbound and 
loosely correlated with m. The calorimetric excitation en- 
ergy sorting leads to comparable results to the multiplicity 
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Table 2. Maximum 72 and Om. values measured in the break known behavior at relativistic energy[5] appears confirmed 

up of a Au system within different data sets sorted in Ztound, in the case of the Nimrod experiment [37| which was per- 

total multiplicity (m) or calorimetric excitation energy (e*). formed with a beam energy as low as 47 MeV/A. This 

Different values for the same case denote different bombarding suggests that a quasi-projectile emission source can be 

energies. Values taken from refs. [5l[lQll3l[33l|34||35||36] extracted^T] in spite of the important midrapidity con- 

iiif=i ii ii — tiibyiiun in the Fermi energy regime [19]. 

72 [5] 72[32]|33]|34] 72 [M] 72 [36] a^M ^^,^6^^ table [3] we can also see that Sl decrease mono- 



Zbound 1-4 1.3 

m 1.85 
1.85 
2.5 



3.2 



3.7 



2.23 



2.5 .12 



Table 3. Maximum 72, and a^n values measured within 
different data sets for various system sizes Zq. Different values 
for the same case denote different targets. Values taken from 
refs.[3M3ZlES] 



< Zo> 


72[35l|37][36] 




o-m[35l|37||36] 


76 


2.5 


1.49 


0.14 


59 


3.7 


0.85 


0.12 


43 


2.4 


0.73 


0.13 


27 


1.75 


0.39 


0.125 


16 


1.19 


0.22 


0.114 




1.17 


0.22 


0.114 




1.16 


0.22 


0.114 



sorting. The value of 72 is slightly increased, which may 
be explained by a reduced correlation of e* with respect 
to m with the total fragment charge, since the excita- 
tion energy contains the extra information of the kinetic 
energy of the fragments. However the effect goes in the 
opposite direction as the fluctuation of is concerned. 
A detailed study of the correlation coefficient between the 
considered observables and the sorting variables is needed 
to fully understand thess trends. It is also possible that 
the fluctuations obtained with these two sortings may be 
compatible within error bars, which stresses the impor- 
tance of an analysis of errors. 

The fluctuation values appear to be largely indepen- 
dent of the reaction mechanism and incident energy [5115^ 
[34] . The only exception is the value 72 « 2.5 obtained from 
emulsion data in ref. |33j , which is significantly higher than 
the values obtained at the other bombarding energies for 
the same system. Such anomaly might be due to the pres- 
ence of fission events that have been excluded in the other 
analyses [32ll34j . The independence on the incident energy 
tends to show that the fragmentation process is essentially 
statistical. 



3.2 Effect of the system size 

The effect of the system size is further analyzed in table 
[3] All presented data are sorted in bins of calorimetric 
excitation energy. 

The fluctuation properties of quasi-projectile decay ap- 
pear to be largely independent of the target. This well 



tonically with the system mass. The evolution with the 
sJ'Stem size, at least in the size range analyzed, appears 
as a simple scaling behavior as shown by the fact that the 
normalization to the source size in am makes the fluctu- 
ation almost independent of the size. Similar conclusions 
can be drawn concerning the 72 observable, even if the be- 
havior for the heaviest sources is less clear. This interesting 
scaling behavior should be confirmed using hyperscaling 
techniques [TTj. 

To conclude, we have seen that fluctuations can vary of 
a factor two changing the sorting variable. This stresses 
the need of confronting the experimental data with sta- 
tistical predictions containing the same constraints, i.e. 
performed in the adapted statistical ensemble. Interesting 
enough, when the same sorting is adopted the different 
available data sets agree within « 15% both in the value 
of the peak and in the position where the peak is observed. 
More data are needed to confirm these trends. 



3.3 Configurational energy fluctuations 

One of the most interesting aspects of studying fluctuation 
observables, is their possible connection with a susceptibil- 
ity or a heat capacity via eq.Q. Conflgurational energy 
fluctuations have been studied at length by the Multics 
collaboration [38. 39. ,40j and by the Indra collaboration[5Ul 
|4T1|421|43] on Au sources. The observable used in these 
studies is an estimation of the energy stored in the con- 
flgurational degrees of freedom at the time of fragment 
formation, deflned as follows 

i=l 

+ Q{Z,,Ai)Mf (10) 

+ V,oui {{ZI],Vfo) 

where Q indicates the mass defects and Vcoui the coulomb 
energy. The measured fragment charges Zi and Icp multi- 
plicities Mi are corrected in each event to approximately 
account for secondary decay 



Zi 



Zf = Z, + {Mf," + 2MZ)- j,^^^ 

l^i=i Zi 

Mf = M, - (M™) 



(11) 
(12) 



where {Mf") is the estimated multiplicity of secondary 
emitted light charged particles for each calorimetric exci- 
tation energy bin. 

Three quantities need to be estimated in each excita- 
tion energy bin to compute Ei 
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Fig. 5. Left part: percentage of secondarily emitted light 
charged particles taken from correlation functions measure- 
ments (see ref.[47)). Right part: total measured fragment ki- 
netic energy (points) compared with Coulomb trajectories cal- 
culations where the volume is changed from 6Vb to 3Vb. Both 
quantities are plotted as a function of the calorimetric excita- 
tion energy. Figure is taken from [39]. 

1. The freeze out volume Vpo which determines the total 
Coulomb energy. Its average value is deduced from the 
measured fragment kinetic energies through Coulomb 
trajectories calculations (see Figure EJ right part). 

2. the average multiplicities of secondarily emitted par- 
ticles (M/^p to account for side feeding effects. They 
are deduced from fragment-particle correlation func- 
tions (see Figure [5l left part). 

3. the isotopic content A^/Zf of primary fragments. It 
is assumed that it is equal to the isotopic content of 
the fragmenting system. This quantity allows in turn 
to determine the number of free neutrons at freeze out 
from baryon number conservation. 

A general protocol has been proposed to minimize the 
spurious fluctuations due to the implementation of this 
missing information[39]. The resulting fluctuation of Ei 
aj — of. is shown for different Multics data in flgure [H 
The temperature has been estimated alternatively using 
isotopic thermometers or solving the kinetic equation of 
state and comes out to be in good agreement [40 with the 
general temperature systematics [44] (around 4.5 MeV in 
the fragmentation region). Similar to the other fluctua- 
tion observables, configurational energy fluctuations show 
a well pronounced peak at an excitation energy around 5 
A.MeV. This general feature is apparent in Multics [40 , 
Indra[43]. Isis[45] and Nimrod[37] data. The only excep- 
tion is EOS datajlg] where this fluctuation appears mono- 
tonically decreasing. 

In the hypothesis of thermal equilibrium at the freeze 
out configuration this fluctuation is a measurement of the 
heat capacity according to eq. ([9|) . The value expected for 
this fluctuation in the canonical ensemble can be written 
as cr^^„ = CfcT^. The kinetic heat capacity Ck is calculated 
from the measured fragment yields [39] . We can see that 
the fluctuation peak overcomes the upper classical limit 
Ck = 3/2 suggesting a negative heat capacity as expected 
in a first order phase transition analyzed in the micro- 
canonical ensemble |481l49] . 

The same analysis performed on Indra data of central 
Xe-|-Sn collisions at different bombarding energies leads 




2.5 5 7.5 2,5 5 7.5 



P/Ao (A. MeV) 

Fig. 6. Left side: normalized fluctuation of Ej and estimated 
Ck (see text) as a function of the calorimetric excitation energy. 
Grey zone; peripheral 35A.MeV Au-|-Au collisions. Symbols: 
central Au+C, Au-t-Cu, Au-|-Au at 25 and 35 A.MeV. Right 
side: heat capacity from eqU) Figure is taken from [40j . 

to compatible temperatures and volumes and a fiuctua- 
tion estimation that agrees within 25% with the presented 
Multics results |39j, as shown for the 32 A.MeV data in fig- 
ure [7| (upper part). In the absence of isotopic resolution 
for fragments, Coulomb repulsion cannot be distinguished 
from a radial collective expansion due to a possible ini- 
tial compression. If an important radial flow component 
is assumed for these central collisions, data can also be 
compatible with a bigger freeze out volume (lower part of 
the figure) leading to a shift of the abnormal fluctuation 
behavior towards lower energy. This volume/flow ambi- 
guity in central collisions can only be solved with third 
generation multidetectors [57] . 

Indra data on the same Au quasi-projectile analyzed 
by the Multics collaboration lead to a fluctuation mea- 
surement about 40% lower, see figure IH This difference 
is tentatively explained as an effect of emission from the 
neck which leads to a reduced occupation of the available 
phase space [13]. 

Recent Nimrod data[33 on the fragmentation of a much 
lighter system show a similar value for the energy corre- 
sponding to the fluctuation peak, but fluctuation a factor 
10 lower than for Multics data, as shown in figure (5] If 
we consider the global fluctuation < A^) > a'f. without 
the normalization to the estimated temperature, this fac- 
tor is reduced to about a factor 4. These results go in 
the same direction as the general behavior of 17^ that we 
have analyzed in section [221 Recall that the fluctuation of 
the biggest fragment for the quasi- Au source [35] is a fac- 
tor 6.8 higher than for the quasi- Ar[37[. This fluctuation 
reduction seems then to be a general feature of light sys- 
tems fragmentation and has been tentatively explained as 
an effect of the higher temperature that light systems can 
sustain[37j. In this interpretation, a higher temperature 
region of the phase diagram, possibly above the critical 
point, is explored in the fragmentation of light systems, 
and the first order phase transition observed in heavy nu- 
clei becomes a smooth cross over. 

As a general remark, the configurational energy fiuc- 
tuation signal is a very interesting one due to its possi- 
ble connection with a heat capacity, but it is also a very 



F. Gulminelli, M.D'Agostino: Fluctuations of fragment observables 



7 



0-2.5 
1 2 



b- 1 



M M 1 M M 1 M M 1 M m 
■ ■ ■ « ■ n 


J M 1 1 1 M 1 1 1 M M 1 1 M 

1 r - 


1 1 i_i 1 i-i-rri 1 1 1 1 1 1 1 1 1 1 





234-56789 10 
THERMAL EXCITATION ENERGY (A.MeV) 




234-56789 10 
THERMAL EXCITATION ENERGY (A.MeV) 

Fig. 7. Normalized fluctuation and kinetic heat capacity 
(stars) for 32 A.MeV Xe+Sn central collisions measured by 
the INDRA collaboration as a function of the calorimetric ex- 
citation energy with two different hypothesis on the freeze out 
volume. The histogram gives the event distribution. Figure is 
taken from [4T1I39] . 



b 



b 



Fig. 8. Normalized fluctuation and kinetic heat capacity 
(stars) for 80 A.MeV Au+Au peripheral collisions measured 
by the INDRA and ALADIN collaboration as a function of the 
calorimetric excitation energy, for all quasi-projectile events 
(left side) and after subtraction of events elongated along the 
beam axis (right side). Figure is taken from |43) . 



indirect and fragile experimental signal which needs pre- 
cise calorimetric measurements, a careful data analysis, 
extensive simulations to assess the effect of the different 
hypotheses in the event sorting and reconstruction pro- 
cedure. Moreover the different techniques to exclude or 
minimize preequilibrium and neck emission seem to have 
a strong influence in the absolute value of fluctuations. 

The evaluation of systematic errors in fluctuation mea- 
surements is necessary to achieve a quantitative estima- 
tion of fluctuations: some first encouraging results in this 
direction have been presented in ref . [40j . The confirma- 
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Fig. 9. Normalized fluctuation and kinetic heat capacity for 47 
A.MeV Argon quasi-projectiles on different targets measured 
by the NIMROD collaboration as a function of the calorimetric 
excitation energy. Figure is taken from [37] . 



tion (or inflrmation) of the fluctuation enhancement is 
certainly one of the most important challenges of the field 
in the next years with third generations multidetectors. 



4 Open questions 

The possibility of accessing a thermodynamic information 
on the nuclear phase diagram from measured fragment 
properties entirely relies on the representation of the sys- 
tem at the freeze out stage as an ideal gas of fragments p6] 
in thermal equilibrium. This is true for fluctuation observ- 
ables as well as for all other thermodynamic analyses [HI 
. This is an important conceptual point which is presently 
largely debated in the heavy ion community. 

A first open question concerns the structure of the sys- 
tems at the freeze out stage, i.e. at the time when frag- 
ments decouple from each other. Contrary to the ultra- 
rclativistic regime [5^, we do not expect much difference 
between the chemical and kinetic decoupling times due 
to the small collective motions implied in these low en- 
ergy collisions. We can therefore speak at least in a first 
approximation of a single freeze out time. If at this time 
the system is still relatively dense, the cluster properties 
may be very different from the ones asymptotically mea- 
sured, and the question arises[T3 whether the energetic 
information measured on ground state properties can be 
taken backward in time up to the freeze out. Calcula- 
tions from classical molecular dynamics [5T] show that the 
ground state Q-value is a very bad approximation of the 
interaction energy of Hill clusters in dense systems. This 
is due both to the deformation of clusters when recognized 
in a dense medium through the Hill algorithm, and to the 
interaction energy among clusters in dense configurations 
where clusters surfaces touch. As a consequence, compara- 
ble fluctuations are obtained in the subcritical and super- 
critical region of the Lennard-Jones phase diagram. This 
result is shown in flg llOl Calculations in a similar model, 
the Lattice Gas model, show that even in the supercritical 
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Fig. 10. Phase diagram of 64 Lennard Jones particles confined 
in a box. Filled triangles give the coexistence border. Isocon- 
tours give values of normalized fluctuations (j'i/<^'can calculated 
from the ground state Q value of clusters defined with the Hill 
algorithm. Figure is taken from [51) . 



regime the correct fluctuation behavior can be obtained if 
both the total energy and the interaction energy are con- 
sistently estimated with the same approximate algorithm 
as it is done in the experimental data analysis [52]. Indeed 
the high value of the estimated configurational energy Q 
fluctuations is essentially due to the spurious fluctuation 
of the total energy Ek + Ei obtained when Ej is esti- 
mated through Q\ such an effect is eliminated if data are 
analyzed in bins of Ek + Q- This calculation is shown in 

figCH 




-10 1-10 1 

Energy/A (e) Energy/A (e) 

Fig. 11. Normalized fluctuations crf./T'^ as a function of energy 
for a system of 216 Lattice Gas particles in the isobar ensemble 
at different pressures. Full lines: exact results. Symbols: estima- 
tion from the ground state Q value of Coniglio-Klein clusters. 
Dashed lines: as the symbols, but data are sampled in bins of 
energy reconstructed from cluster kinetic energies and sizes. Ac 
gives the critical pressure. Figure is taken from |52) . 



A second related question which needs further work 
is the relevance of the equilibrium assumption at freeze 
out. Molecular dynamics models applied to study the time 
evolution of the reaction[2T ll531l54p 55j predict that the de- 
coupling between fragment degrees of freedom (freeze out) 
occurs very rapidly during the reaction. At this stage how- 
ever the configuration is considerably diluted due to the 
early presence of collective motions pT[|20] . An example 
taken from classical molecular dynamics for an initially 
equilibrated compact configuration freely evolving in the 
vacuum, is shown in fig |12l At this reaction stage cluster 
energies may be well approximated (within a side feeding 
correction) by their asymptotically measured values, but 
it is not clear whether this configuration can correspond 
to an equilibrium, more precisely whether the hypothesis 
of equiprobability of the different charge partitions holds. 




100 200 3 00 100 200 300 30 100 130200 230 300 

tiLnefTijJ time (i^) time (t^) 

Fig. 12. Time evolution of a Lennard Jones system initially 
confined in a dense supercritical configuration and freely ex- 
panding in the vacuum at different total energies. Upper part: 
MST fragment size distribution at different times. Lower part: 
average kinetic energy (full lines), total (lower symbols) and 
normalized (upper symbols) kinetic energy fluctuations, and 
size of the largest MST cluster (dashed lines). Abnormal fluc- 
tuations in these units correspond to Ak >~ 0.7. Figure is 
taken from 1531. 



5 Conclusions and outlooks 

In this paper we have presented a short review of the ex- 
perimental as well as theoretical studies of fluctuation ob- 
servables of fragments produced in a multifragmentation 
heavy ion reaction. The aim of these studies is the under- 
standing of the nature of the nuclear fragmentation transi- 
tion as well as the thermodynamic characterization of the 
finite temperature nuclear phase diagram. This vast and 
ambitious program is still in its infancy. Many promising 
results already exist, but the analyses are not yet conclu- 
sive and need to be intensively pursued in the future. 
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The nuclear fragmentation phenomenon, well docu- 
mented by a series of independent experiments [7 , presents 
many features compatible with a critical phenomenon [1 1 j 
or a phase transition [9, 25]. Only a careful study of fluc- 
tuation properties will allow to discriminate between the 
different scenarii. Even more important, the phase dia- 
gram of finite nuclei is theoretically expected to present an 
anomalous thermodynamics [49,48^ which should be char- 
acteristic of any non extensive system undergoing first or- 
der phase transitions in the thermodynamic limit. Once 
the difficulties linked to the imperfect detection and sort- 
ing ambiguities will be overcome, fluctuation observables 
will be a unique tool to quantitatively study this new ther- 
modynamics with its interdisciplinary applications [49il48l 
[56]. 

From the theoretical point of view, the theoretical con- 
nections between fluctuations and susceptibilities in the 
different statistical ensembles are well established, and 
the different experimental contraints can be consistently 
adressed by the theory. However, the evaluation of a ther- 
modynamics for a clusterized system opens the difficult 
theoretical problem of cluster definition in dense quan- 
tum media. To produce quantitative estimations of mea- 
surable fluctuation observables the pertinence of classical 
models has to be checked through detailed comparisons 
with microscopic [54] and macroscopic [26] nuclear models. 

On the experimental side, multiplicities and size fluctu- 
ations agree reasonably well if comparable size fragment- 
ing systems are studied, even if the effect of the system 
size has to be clarified. Configurational energy fluctuations 
are especially interesting because of their possible connec- 
tion with a heat capacity measurement. The methodology 
to extract such fluctuations from fragmentation data is 
presently under debate, in particular a careful analysis 
of systematic errors is presently undertaken[40]. From a 
more conceptual point of view, the influence of the differ- 
ent time scales in the reaction dynamics has to be clari- 
fied. Configurational energy fiuctuations may be subject 
to strong ambiguities since they use information from all 
the particles of the event, and this information is inte- 
grated over the whole reaction dynamics. In this respect, 
an interesting complementary observable may be given by 
fluctuation of the heaviest cluster size [25, 29:. 

To solve the existing ambiguities we need full compar- 
isons with a well defined protocol and consistency checks 
between different data sets. The simultaneous measure- 
ment of fragments mass and charge on a 47r geometry [57] 
will be essential to measure the basic variable of any ther- 
modynamic study, namely the deposited energy. No defini- 
tive conclusion about the occurrence of a thermodynamic 
phase transition and its order can be drawn without this 
detection upgrade. 
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